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ABSTRACT 



Electron accelerations at high Mach number collision-less shocks are investi- 
gated by means of two-dimensional electromagnetic Particle-in-Cell simulations 
with various Alfven Mach numbers, ion-to-electron mass ratios, and the upstream 
electron j3e (the ratio of the thermal pressure to the magnetic pressure). We 
found electrons are effectively accelerated at a super-high Mach number shock 
[Ma ~ 30) with a mass ratio of M/m — 100 and — 0.5. The electron shock 
surfing acceleration is an effective mechanism for accelerating the particles to- 
ward the relativistic regime even in two dimensions with the large mass ratio. 
Buneman instability excited at the leading edge of the foot in the super-high 
Mach number shock results in a coherent electrostatic potential structure. While 
multi-dimensionality allows the electrons to escape from the trapping region, they 
can interact with the strong electrostatic field several times. Simulation runs in 
various parameter regimes indicate that the electron shock surfing acceleration 
is an effective mechanism for producing relativistic particles in extremely-high 
Mach number shocks in supernova remnants, provided that the upstream elec- 
tron temperature is reasonably low. 

Subject headings: acceleration of particles - cosmic rays - plasmas - shock waves 
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Introduction 



Origin of cosmic rays has been a longstanding issue in astrophysics. A rapid progress 
in understanding their accelerators has been made by observations of radio and X-ra y 



synchrotron emissions from the shell of supern o va remnants (SN Rs) ( iKoyama et al 



1995 



Winkler fc Long 



1997 



Bamba et al. 



2003 



Long et al.l 



20031 ). Those observations 



have indicated that electrons are strongly accelerated up to TeV energies out of thermal 
population in the vicinity of non-relativistic collision-less shocks. 

Diffusive shock acceleration (DSA) is a widely accepted theory responsible for 
particle acceleration at sh ocks, which nicely predicts a power low energy spectrum (e.g., 



Blandford &: Eichler 



19871 ) . This theory has a chieved many successes in explaining energized 



ions, for example, at the earth's bow shock ( IScholer et al. 



19801 ). However, it cannot be 



straightforwardly applied to electron accelerations at SNRs: In DSA process, particles are 
scattered by MHD waves back and forth across the shock. The process is thus unlikely to 
occur for thermal electrons whose gyro radius is much smaller than the wavelength of the 
MHD waves or the shock width (ion inertia length). They are strongly tied to the magnetic 
field line and convect toward downstream adiabatically. This small-scale property of the 
electrons makes it challenging to explain the electron accelerations by DSA at SNR shocks. 
A pre-acceleration mechanism for the electrons that enables DSA to work subsequently has 
been proposed for the last decade in order to overcome this difficulty, which is now known 
as the injection problem. 

In supercritical quasi-perpendicular shocks (M4 > 2.7, where is the Alfven Mach 



number), it has been p ointed out that plasma microscopic instabilities play 



in energizing electrons ( 



2003 



Papadopoulos 



Scholer fc Matsukiyc 



1988 



Cargill fc Papadopoulos 



1988 



important ro. 



Scholer et al. 



es 



20041 ). In the particular case of high Mach number shocks. 



electrostatic instabilities arise due to velocity difference between the reflected ions and 
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the i ncident (upstream) ions and electrons (jPapadopoulod Il988l : ICargill &: Papadopoulos 



19881 ). Buneman instability is initially excited by the reflected ion beam interacting with 
the upstream electrons. The instability strongly heats the electrons perpendicular to the 
magnetic field, which then gives a favorable condition for destabilizing the ion acoustic 
instability in the foot region. Particle-in-Cell (PIC) simulations, which follow ion and 
electron motions along with the electromagnetic fields development in a self-consistent 
manner, have confirmed those ideas. Furthermore, a rapid acceleration of the electrons 



instability was found to contribute to nonthermal e 



2000 



McClements et al. 



2001 



Hoshino fc Shimada 



2002 



production ( 


SI 


limada & Hoshino 


Shimada et al. 


20101). This 



electron shock surfing acceleration has been considered to be a very efficient mechanism at 
extremely-high Mach number shocks such as SNR shocks that injects thermal electrons to 
DSA. 

It is only recently that evolutions of no n-relativistic shock s has been investigated by 



means of two-dimensional PIC simu 



Lembege et al. 



2009 



Kato fc Takabe 



201G 



fiUmeda et al. 


2009; 


An 


lano & Hoshino 


2009a 




Riquelme & Spitkovskv 


2011 


). The efficiency of 



the electron shock surfing acceleration in multi dimensions has been discussed in different 
physical parameter regimes, such as the ion-to-electron mass ratio, Alfven Mach number, 
upstream electron /3e (the ratio of the thermal pressure to the magnetic pressure), and 
magnetic field configurations; the two-dimensional PIC simulation still requires large 
computational resources and only allows limited simulation runs. Consequently, the 
results lead to different conclusions. The simulation study with a relatively low mass ratio 
(M/m = 25) concluded that the shock surfing acceleration is an effec tive mechanism eve n 



in two d imensions in a h i gh (M a = 14) Alfven Mach number shock (lAmano &: Hoshino 



2009ah . 



Kato &: Takabd ( 120101 ) examined an extremely-high Alfven Mach number 



[Ma = 130) shock with a high upstream electron /3e condition. In their simulation 
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result, the downstream electron energy spectrum exhibited almost a Maxwellian. They 
attributed the lack of the nonthermal population to the in-plane mag netic field configurat ion 



which w a s different from the out - of-pla ne configuration adopted by 



Amano fc Hoshino 



( l2009al ). iRiquelme fc SpitkovskyI ( 1201 ll ) conducted simulation runs with various mass 
ratios and concluded that the Buneman instability is suppressed in a large mass ratio case 
(M/m = 400) under the condition of Ma = 7 and the electron shock surfing acceleration 
will not be realized at SNR shocks. Those contrary conclusions regarding the electron shock 
surfing acceleration at SNRs are mainly due to the limited choices of the mass ratio, Alfven 
Mach number, and the upstream electron /3p , which should be carefully chosen based on the 



theoretical estimations ( jPapadopoulos 



1988 



Matsukiyo 



2010|). 



In this paper, we present two-dimensional electromagnetic PIC simulations of collision- 
less shocks with various mass ratios, Alfven Mach numbers, and the upstream electron /3e 
conditions in order to clarify the efficiency of the electron shock surfing acceleration at SNR 
shocks. The parameters are carefully chosen in accordance with the linear and nonlinear 
theories of the Buneman instability in the foot region, which are reviewed in Section 2. 
In Section 3, the numerical setup as well as our two-dimensional PIC code are described. 
In Section 4, we first overview a simulation result in which the maximum energy gain of 
the electrons is observed. Electron acceleration processes are analyzed in detail. Then 
dependence of the acceleration mechanisms on the physical parameters are discussed. The 
results are summarized and discussed in the last section. 



2. Unstable condition and the saturation level of Buneman instability in the 

foot region 

The electron shock surfing acceleration is basically supported by the electrostatic field 
excited by the Buneman instability at the leading edge of the foot region. The free energy of 
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the instability is the relative drift motion between the incoming electrons and the reflected 
ions. The region becomes unstable when the relative speed exceeds the thermal speed of the 
electron, ^V/vthe > 1- In applying the unstable condition to the foot region, the velocity of 
the incoming ion and electrons, and the reflected ions are expressed as 



Vr = -Vo, 

e — :r- — i^o, 
1 + a 



(1) 
(2) 
(3) 



by assuming the zero net current condition 

rioVi + TirVr- {no + nr)Ve = 0, (4) 

where no is the upstream number density, rir is the number density of the reflected ion, Vq 
is the upstream speed, Vi is incoming ion speed, Ve is the incoming electron speed, and Vr 
is the reflected ion speed, a = n^/no is the density ratio of the reflected and incoming ion s, 



which is typically ^ 0.2 for high Mach number shocks (jPapadopoulos 



1981 



1988 



1982 



Leroy 



Leroy et al. 



19831 ). Using Al^ = Ve — Vr, the unstable condition reads as 



1 + a ^fM\ 
2 \m J 



(5) 



where Ma is the Alfven Mach number, M and m are the ion and electron rest mass, and /9e 
is the ratio of the electron pressure to the magnetic pressure (plasma beta). 

Given the Buneman instability being excited, the electrons are accelerated by the 
convective electric field while they are trapped in the electrostatic potential well. The 
maximum energy gain is determined by force balance between the trapping force by the 
saturated Buneman instability {qEBi) and the Lorentz force {qvBo/c) for escaping. To be 
accelerated to ward relativistic speeds (v ~ c), the electric force must be larger than the 



Lorentz force (IHoshino fc Shimada 



2002): 



E 



BI 



> 1. 



(6) 
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The saturation level of the Buneman instability {\Ebi\) can be estimated by introducing an 
energy conversion rate C of the electron drift energy density as 

^ = l-mnoAV'C. (7) 

'he conversion rate C has been analyzed by the nonlinear theory of the Buneman instability 



( llshihara et al. 



19801 ) and found to depend weakly on the mass ratio C ~ (m/M)^/'^. 
In multi dimensions, the saturation level is reduced to ~ 25% of above estimate due 
to the resonant waye-part icle interactions with waves oblique to the beam direction 



( lAmano fc Hoshino 



2009bl ). Using AV = Ve — Vr and eq.(I7]), and taking account of the 



multidimensionality, eq. (|6]) leads to 

/ M\ ^ 

MA>{l + a)i — ) . (8) 

Eq. ([5]) indicates that the electron shock surfing acceleration occurs at a sufficiently 
high Alfven Mach number shock in a low beta plasma. Furthermore, for non-thermal 
relativistic electrons to be observed, eq. (|8]) shows that the mass ratio restricts Alfven 
Mach number more severely with a power of 2/3 rather than 1/2 in the unstable condition, 
indicating a need for numerical simulations of super-high Alfven Mach number shocks with 
larger mass ratios, for instance, > 26 for M/m = 100. In other words, simulations with 
artificial small mass ratios would easily give an unrealistic electron acceleration in relatively 
low Alfven Mach number shocks. 



3. Numerical setup 

Two-dimensional electromagnetic Particle-in-Cell simulation code is used for examining 
shock evolutions. The code is parallelized via domain decomposition by using Message 
Passing Interface (MPI) library and OpenMP and runs efficiently on massively parallel 
supercomputer systems. 
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Particles are continuously injected from one boundary (x = 0) with a speed Vq toward 
the other end {x = L^) where the particles are reflected (injection method), resulting in 
the shock formation that propagates toward x = 0. Number density in the upstream is 
40 particles per cell for each species (ion and electron). The injected plasma carries the z 
component of the magnetic fleld (i?o) and the convective electric fleld Ey = Eq = Vq-Bq- 
Thus we deal with purely perpendicular shocks in the downstream-rest frame. The periodic 
boundary condition is applied in the y direction. Simulation box sizes in the x and y 
directions are = 12 Vo/Qgi = 12 M^Aj and Ly = 5 A,, respectively, where Ma = Vq/Va, 



^Igi = qEq/Mc is the ion gyro frequency, Aj = c/upi = m j ^-nns^q^ is the ion inertia length, 
Va = Bq/ y/AirMno is the Alfven speed, and riQ is the number density in the upstream region 
(40 particles per cell). Other notations of c and q respectively represent the speed of light 
and elementary charge. The grid size Ah and the time step size At are set as Ah = Xd 
and ujpeAt = 0.025, where Xd = a/T^T^^o^ is Debye length in the upstream region, 
is the electron temperature, and oope = ^yArcnQq"^ /m is the electron plasma frequency. In 
the following, we use the units of c/ojpi, no, Bq, c, and mc^ for time, space, number 
density, electromagnetic flelds, velocity, and energy, respectively, unless otherwise stated. 

We have carried out several runs with various ion to electron mass ratios (M/m), the 
Alfven Mach numbers (Ma = Vq/Va), and the electron plasma f3e = SirnoTe/ Bq, while the 
ratio of the electron plasma to gyro frequencies (upe/^ge) and the temperature ratio of the 
ion to the electron [Ti/Tg) are flxed to Upe/^lge = 10.0 and Ti/T^ = 1, respectively. Two 
simulation runs use the mass ratios of 25 and 100 for flxed values of Ma = 10 (respectively 
Ma ~ 14 and 16 in the shock-rest frame) and /3e = 0.5 in order to examine the mass ratio 
dependence of th e electron accelerati o n (Run s A and B). Parameters in Run A are the same 



as used those by 



Amano &: Hoshind (l2009al ). Another simulation run uses M/m = 100 



with Ma = 20 {Ma ~ 30) for comparison with the lower Alfven number cases (Run C). 
We also examined a hot electron case with = 4.5 for comparison with the high Alfven 
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Mach number and a lower electron temperature case (Run D). We totally examined four 
simulation runs (Runs A-D) with various M^-, M/m, and /3e to assess importance of these 
parameters to electron acceleration mechanisms. The upstream parameters of the four 
simulation runs are summarized in Table [1] The largest computational resources are used 
in Run C, in which 5 x 10^ particles are followed in a simulation domain with 24001 x 1024 
grid points. 



4. Simulation Results 



4.1. Electron acceleration in Run C 



4.1.1. Overview of 2D shock structure 



We first present an overview of a two-dimensional perpendicular shock structure in 
Run C in which the maximum energy gain of the electron was observed. In the following 
results, the data at a time of figjT = 6 is used, after which the energy spectrum of the 
electron in the downstream region is not significantly changed. 

Figure [Il^a) shows a two-dimensional spatial profile of the z component of the magnetic 
field. The structure is characterized by the foot region which extends from X = 63 c/upi 
to the ramp and the overshoot around X = 74 c/upi, and the downstream region in 



X > 74 c/ojpi. The width of the foot (the sh ock transition region) is we. 



radius of the reflected ion (L ~ 0.7 Vo/VLgi) ( iLeroy 



1983 



Shimada et al. 



scale d by the gyro 



2010( 1. 



Figure [^b) shows the phase space plot for ions. The ions in the foot consist of the 
incoming cold population and the relatively hot ion moving toward the upstream region 
which is specularly reflected by the shock potential at X ~ 74 c/upi. At the leading edge 
of the reflected ion, electrons are preferentially accelerated immediately after entering the 
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foot as seen in the energy-space diagram in Figure [U^c). Just behind the region around at 
X ~ 63 c/ujpi, the maximum energy of the electrons increased to 7 ~ 2. Then they are 
further heated in X > 66.5 c/upi in the foot. In the vicinity of the overshoot, the electron 
energy again increased up to 7 ~ 9. 

Zoomed-up spatial profiles surrounded by the white dashed line in Figure [U^a) are 
shown in Figures [T](d) and (e). At the edge of the foot, we can identify fine scale structures 
of the electron number density in Figure [T](d), whose spatial size is much smaller than 
the ion inertia scale, i.e., the electron scale that is embedded in the MHD scale shock 
structure. This fine structure is accompanied by electrostatic fields. Figure [^e) shows 
the electrostatic field strength |-Eest| = | — V0| given by solving the Poisson equation of 

= — V ■ E. The amplitude reaches twice as large as the background magnetic field 
(-Bo), which is much larger than the upstream convective electric field \Ey/ Bq\ = Vq/c = 0.2. 
Coherent wave trains lie in the x-y p lane rather than exhibiting two-dimensional isolated 



structures flAmano fc Hoshino 



|2009aJ). This strong electrostatic field at the edge of the foot 
is caused by an linear instability between the incoming electron and the reflected ion, which 
one usually refers to as Buneman instability. The Buneman instability heats the electrons 
perpendicular to the magnetic field and sets a favorable condition for the ion-acoustic 
instability, which further energized the electrons in X > 66.5 c/upi in the foot. 

In order to identify the instability in the two-dimensional plane as the Buneman 
instability, we focus on the closed up region of interest in Figures [U^d) and (e). Figure [2](a) 
shows the zoomed up profile of the electrostatic potential energy q(f). In this region, the 
cold upstream and the reflected gyrating ions constitute the distribution function shown 
in Figure [2]^b). The number of reflected ions with respect to the upstream po pulation is 



nr/rip = C( = 0.22 as expected in previous high Mach number shock simulations (ILeroy et al. 
IQSalLerovlllQSaf ). 
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The relative drift speed between the two populations is AV ~ 0.4c. This velocity 
difference is much larger than the upstream electron thermal speed of 0.07c, thus satisfying 
the unstable condition of the Buneman instability (eq. Figure |2|^c) shows the Fourier 

amplitude of the electrostatic potential of Figure [2]^a) in the k^-ky space in which Hann 
window is applied in the x direction. One can identify a strong peak appearing around 
{kx,ky) = (0.6,0.3) Upe/AV. This peak corresponds to the coherent wavy structure that 
obliquely lies in the x-y plane and \k\ = ^/k^ + k^ = 0.7 Upe/AV agrees with the fastest 
growing mode at k^AV/Qp f. ~ 1 predicted by the linear analysis of the Buneman instability 



(lAmano fc Hoshino 



2009bh . 



4.1.2. Electron acceleration mechanisms 

In Run C, a high energy tail of the electron energy distribution function which extends 
to 7 ~ 9 is observed. Figure [3] shows the energy spectra of the electron in the foot and 
downstream regions. Particles are sampled in regions 64 c/upi < X < 71 c/upi for the foot 
spectrum and 74 c/upi < X < 88 c/upi for the downstream spectrum. In the foot region, 
non-thermal electrons are already produced. The energy distribution consists of the cold, 
non-accelerated and accelerated populations. The highest energy in this region is 7 ~ 3.5. 
In the downstream region, a clear two-component distribution is found. One is the thermal 
electrons whose temperature is Te/mc^ = 0.13 and the other consists of accelerated particles 
whose energy reach 7 ~ 9. Acceleration mechanisms responsible for these non-thermal 
electrons are found in two ways which are described in the following. 

First, the electrons are accelerated at the edge of the foot by the shock surfing 
acceleration as manifested in Figures [U and El When a particle entered the unstable region 
of the Buneman instability, it is accelerated by the convective electric field while being 
captured by the electrostatic potential. The present two-dimensional simulation shows 
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in Figure Hl^a) that a sampled electron can be escaped from the potential well at a time 
after gaining a fraction of energy and drifts toward the downstream. Since the length 
of the potential well is limited to the order of the ion inertia length, the particle can 
escape from the well even though the electric field strength is sufficiently large. This is 
the multi-dimensionality effect. Nevertheless it can enter the unstable region from the 
downstream-side again since the gyro radius of the accelerated particle becomes larger 
than before. Then the second surfing acceleration is realized. In the trajectory in Figure 
Ht^b), the energy gain is more effective than in the previous encounter, and the energy 
reached 7 ~ 2. The third encounter with the unstable region is possible in the same 
manner, and the energy is increased to 7 ~ 3 (Figure Hl^c)). For accelerated particles, they 
typically experience the shock surfing acceleration three times as shown in this example. 
The present acceleration mechanism can be applied to lucky particles that will constitute 
the non-thermal population, and most of the particles that will constitute the thermal 
population drift toward the downstream without being captured as shown with cyan lines 
in Figure IH Nevertheless the shock surfing acceleration is so effective as it can energize the 
electrons from 7 ~ 1 to the relativistic regime of 7 ~ 3 within a time scale of ~ 20 ^ge^ or 
~ 0.2 n~l in this particular case with M/m = 100. 

Time histories of the energy and the magnetic moment (the first adiabatic invariant) 
normalized by Uq/2Bq = Vq^/2(1 — Vq/c'^)Bq are shown in Figured The energy history 
of the accelerated particle is threefold corresponding to Figures |4]^a)-(c). The maximum 
energy gain is found at QgeT ~ 517 when 7 — 1 is rapidly increased from 0.05 to 1.0 within 
a time scale of 5 Qg^ (Figure |5]^a)). Since the shock surfing acceleration is due to the 
direct acceleration by the DC electric field, the magnetic moment also rapidly increases as 
the energy is increased (Figure [St^b)). After the third surfing acceleration at QggT ~ 530, 
the energy and the magnetic moment are kept constant until QgeT = 550. On the other 
hand, the non-accelerated particle remains at the initial energy level while the normalized 
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magnetic moment fluctuates around one; neither obvious non-adiabatic nor adiabatic 
accelerations works on this particle as it travels in the foot region. 

The maximum energy achieved by the surfing acceleration at the edge of the foot is 
up to 7 ~ 3, which is not so sufficiently large as the maximum energy of 7 ~ 9 in the 
downstream (see also Figure [H])- An additional acceleration process is also found as the 
particle approached the shock surface. Figure [6](a) shows a trajectory of the accelerated 
electron in Figure H] at QgeT = 580 after the first acceleration process. The particle 
undergoes the grad-B drift in the -y direction while being deflected at the overshoot. 
Figure MJo) shows that the acceleration starts when the particle encountered with the 
shock surface at flgeT ~ 560. After this time, the energy variation exhibits a sawtooth 
pattern rather than a sinusoidal one in QgeT = 540 — 560; the energy increases rapidly and 
decays gradually. During the acceleration, the magnetic moment stays at a certain level 
on average, but is highly oscillatory (Figure [6](c)). Figure [6](d) shows the running average 
of the magnetic field at the particle's position. The particle experiences larger magnetic 
fields as it is energized in time. These overall signatures in Figures |6]^b)-(d) indicate that 
the acceleration at the shock surface is basically an adiabatic process. 

Figures [Tl^a) and [7](b) shows the magnetic and electric fields' profiles with respect to 
the overshoot position averaged in the y direction (solid lines). The region in the vicinity of 
the overshoot is enlarged. The electric field is calculated in the shock-rest frame based on 
the electron frozen-in condition (the generalized Ohm's law including the hall term) as 

E'yix) = Ey{x) - V,hB,{x), (9) 

where Vsh is the electron convective velocity in the x direction at the overshoot position. 
The amplitude of the magnetic field at the overshoot reached twenty times as large as 
the upstream value within a scale of the ion inertia length, giving rise to the fast grad-B 
drift speed of the relativistic electron. Since the direction of the convective electric 
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field is opposite to the drift motion, tlie electron can gain energy as it travels in the y 



direction. This is similar to the t ypical (ion) shock dri : 



electron ([Armstrong et al. 



1985 



Krauss-Varban et al. 



t acc eleration, but for the relativistic 



19891 ). Furthermore, the electron is 



accelerated by the strong negative Ey in the ramp {X = —0.1 c/upi) unlike the typical shock 



drift acceleration. This electri c field is induced by time yariations of t 



that is, the shock reformation (iLembege fc Savoini 



1992 



Scholer et al. 



le shock structure, 



2003 



Lembege et al. 



20091 ). Figure [7|^c) shows the z component of —V x E = dB/dt. It takes large positive and 



negative values in the ramp in the shock-rest frame, indicating that a new shock front is 
forming ahead of the existing overshoot position. Since the inductive electric field strength 
is larger than the convective electric field (l-Emdl ~ 0.3-Bo), the betatron acceleration is also 



effective for the relativistic electron in the present simulation run. 



4.2. Mass ratio and Alfven Mach number dependencies 

We have shown the shock surfing acceleration and additional adiabatic accelerations at 
the shock surface in Run C in which the maximum energy gain of the electron is observed. 
Here we examine these mechanisms in the other simulation runs with various ion-to-electron 
mass ratios (M/m) and Alfven Mach numbers (M^). 

Figure [8] shows a comparison of the electrostatic field strength in the leading edge of 
the foot. In Run A with M/m = 25 and Mj^ ~ 14, the strong electrostatic field is excited 
indicating the Buneman instability (Fig. |8]^a)). The spatial profile is not like a coherent 
wave form in Run C (Fig. EJ^c)), but consists of isolated peaks distributed in the x-y plane 



( lAmano &: Hoshindl2009al ). The amplitude of the electrostatic field reaches \Eest\/Bo ~ 2.4 
which satisfies the trapping condition of the electron (eq. ([H]))- For a reduced mass ratio, 
the condition of the acceleration toward the relativistic regime (eqs. and ([7])) is relaxed 
to a relatively low Alfven number case, for instance, M/m = 25 requires M^ > 10.3. Run 



- 15 - 



A is indeed the case since Ma ~ 14. 

In Runs A and C, the Buneman instabihty is highly activated and the electrons 
are accelerated in the foot since the trapping condition is satisfied in both simulation 
runs. However, difference arises in number and the maximum energy of the non-thermal 
population in the downstream. We have shown in the previous subsection that the second 
acceleration process is found at the shock surface in Run C which further energizes the 
pre-accelerated electron from 7 ~ 3 up to 7 ~ 7 adiabatically. The key for this acceleration 
is the magnetic and electric fields' profiles around the overshoot. Figure [Tl^a) compares 
the magnetic field structures in Runs A and C. The peak amplitude of the magnetic field 
at the overshoot is Bz/Bq ~ +8.0 in Run A which is much smaller than the peak of 
Bz/Bq ~ +20.0 in Run C. The inductive electric field in the ramp found in Run C is not 
evidently found in Run A (Figures [^b) andl^c)) and thus the betatron acceleration by 
the reformation process is not effective. The pre-accelerated relativistic electron in Run A 
cannot gain much energy through the adiabatic accelerations at the shock surface. 

If the mass ratio is increased with the Alfven Mach number being fixed, one expects the 
Buneman instability to be weakly saturated, or even completely stabilized from eqs. (|5]) and 
([7]). In Run B, M/m = 100 and ~ 16, and therefore the trapping condition > 25.8 
is not satisfied while the unstable condition Ma > 4.2 is realized (Tabled]). Figure [HH^b) 
refiects these conditions that the amplitude of the electrostatic field is \Eest\/Bo ~ 0.5 even 
though the Buneman instability seems to become unstable. Thus the electrons can easily 
escape from the potential well after moderate energization, and the efficient acceleration to 
the relativistic regime is not expected in this run. 
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4.3. Effect of upstream electron temperature 

Upstream electron temperature also affects the activity of the electrostatic field in 
the foot; the Buneman instability is expected to be weakly saturated or even completely 
stabilized in a high electron /3e condition. To assess this effect, we simply increased /3e 
from 0.5 in Run C to 4.5 with M/m and Ma being fixed in Run D. In this case, both the 
trapping condition Ma > 25.8 and the unstable condition Ma > 12.7 are satisfied (Table 
[T]). The resultant electric field profile in the leading edge of the foot is shown in Figure 
[Ht^d). The amplitude of the electrostatic field is so weak {\Eest\/Bo ~ 0.5) that any coherent 
waves are not recognized out of the thermal noise of the hot electron. As in the case of Run 
B, the non-thermal electrons are hardly found in the downstream. 



4.4. Energy spectra in the downstream region 

Lastly, we summarize our simulation runs by showing the energy spectrum of the 
electron in the downstream region. Particles are sampled in a rectangular region of 
0.7 Vo/Qgi X Ly behind the overshoot. Figure [9](a) compares the results from Runs A-C. 
The energy spectrum in Run C exhibits a high energy tail which deviates from the thermal 
population and reaches 7 ~ 9 as was shown in Figure [3l A similar non-thermal electron 
production is observed in Run A in which the Buneman instability is highly activated at 
the edge of the foot as well. These simulation runs satisfy the trapping condition of eq. ([6]). 
When the trapping condition does not hold by increasing the mass ratio, the non-thermal 
population is significantly reduced. This is clearly seen in Run B when compared with 
Run A. It is also notable that fitted gaussian profiles give the electron temperatures of 
Te/mc^ = 0.062 for Run A and Tg/mc^ = 0.033 for Run B. Although the Mach number is 
the same in both runs, the downstream condition is greatly altered by the activity of the 
Buneman instability in the foot which heats the electrons perpendicular to the magnetic 
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field. 

Figure MJo) shows the effect of the upstream electron temperature on the non-thermal 
electron production by comparing the results from Runs C and D. Even though both the 
unstable and trapping conditions are satisfied in Runs C and D, the non-thermal population 
in the downstream in Run D is also significantly reduced. Thus the upstream electron 
temperature (/3e) must be another important parameter that determines the efficiency of 
the electron accelerations. Fitted gaussian profile for Run D gives the electron temperature 
of Te/mc^ = 0.11 which is slightly smaller than in Run C. 

The results are summarized in Table [1] in relation to the upstream conditions. 



5. Summary and Discussion 

We have examined electron acceleration mechanisms at high Mach number shocks by 
means of two-dimensional Particle-in-Cell (PIC) simulations with various ion-to-electron 
mass ratios, Alfven Mach numbers, and the upstream electron Pe- We found electrons are 
effectively accelerated at a super-high Mach number shock (M^ ~ 30) with a mass ratio of 
M/m = 100 and (3^. = 0.5. The electron shock surfing acceleration is found out to be an 
effective mechanism for accelerating the particles toward the relativistic regime even in two 
dimensions with the large mass ratio. Buneman instability in the super-high Mach number 
shock resulted in a coherent electrostatic potential structure which enabled an efficient 
trapping of the electrons. While multi-dimensionality allows the electrons to escape from 
the trapping region, they can interact with the strong electrostatic field several times. This 
multiple interaction e nables effective shock surfi ng accelerations as has been found in the 



small mass ratio case (lAmano fc Hoshino 



2009af ). 



The conditions of the electron shock surfing acceleration toward the relativistic 
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regime have been derived from one- dimensional arguments ( jPapadopoulosI Il988 



Cargill fc Papadopoulos 



1988 



Ishihara et al 



1980l ). Those simple estimations still hold in 



the present two-dimensional simulations as summarized in Figure [TUJ While all our four 
simulation runs satisfies the unstable condition of the Buneman instability (dashed lines), 
the shock surfing acceleration becomes effective in two simulation runs (denoted by "Run 
A" and "Run C" in Figure [TU]) whose upstream conditions are also well above the trapping 
condition of the electron (solid line). A similar aspect holds in recent two-dimensional PIC 
simulations with different parameters from our simulation runs. For example, the electron 
shock surfing acceleration does not play important roles in producing non-thermal electrons 
in shocks with Ma = 7, M/m = 100, 400 an d /3e = 0.5 in which the trapping condition 



does not meet (IRiquelme fc Spitkovsky 



201l|). 



Exception was found for high electron /3g cases. In our high /3e run (denoted by "Run 
D" in Figure [To]), the Buneman instability is destabilized in the foot region. However, its 
peak amplitude is not so large that electrons can be escaped from the trap ping region befo re 



reachi ng the relativistic regime. A similar conclusion was also drawn by .Kato fc Takabe 



fl2010h in which Ma ~ 130, M/m = 30, and (3e = 26 are used (see also Figure [lOD. Although 
their linear analysis revealed that the foot region in their simulation result was unstable to 
the Buneman instability, the resultant energy spectrum in the downstream region showed 
almost a Maxwellian like what we see in our high /3e run. These results indicate that 
we cannot simply understand the high electron /3e simulations from Figure [10] which is 
based on the linear and nonlinear theories of cold plasma, and detailed analysis of the 
saturation mechanism of the Buneman instability with finite electron temperature effects 
is necessary. Extrapolating from the recent and present two-dimensional PIC simulations. 
Figure [10] indicates that the electron shock surfing acceleration is an effective mechanism 
for producing relativistic particles in extremely-high Mach number shocks in supernova 
remnants (SNRs), provided that the upstream electron temperature is reasonably low. 
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In addition to the shock surfing acceleration in the foot, the adiabatic accelerations 
are found to be effective for the pre-accelerated, relativistic electrons in the super-high 
Mach number shock. The process is a combination of the shock drift acceleration and 
the betatron acceleration by the time variation of the ramp due to the shock reformation 
process. Since the mechanism itself is an adiabatic one, the maximum attainable energy by 
this process is limited by the amplitude of the magnetic field at the overshoot: 

^ = (10) 



where the subscripts "/t" and "of " denotes the foot and the overshoot, respectively. 7/t is 
the energy after the shock surfing acceleration in the foot. The magnetic field strength at 
the leading edge of the foot does not strongly depend on the Alfven Mach n umber a nd ou r 



Lerovl fll983h . 



simulation runs exhibit -B/t ~ 2Bq which agrees with the multifiuid model of 
In the Leroy's model, the field strength at the overshoot was also estimated to leading order 

as 

Bov^aBoMl (12) 

where we assume a is a constant. This scaling can be applied to our simulation runs up 
to Ma = 30. For instance. Boy = S.OBq for Ma = UA in Run A and B^v = 20.05o for 
Ma = 30.0 in Run C give a ~ 0.4. Eqs. flTTl) and f[T^ are combined to give the maximum 
energy attained by the adiabatic acceleration around the overshoot as 

^ov = \/0:2Mj7;„ (13) 

which suggests that the electrons accelerated in the foot can be further energized by an 
order of magnitude in extremely-high Mach number shocks in SNRs. Note, however, that 
the high energy particles reside in the vicinity of the shock surface and their energy is 
decreased as they are transmitted in the downstream far away from the overshoot owing to 
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the conservation of the magnetic moment. Non-adiabatic processes, such as a pitch angle 
scattering of the particles by electromagnetic waves at the shock surface, are necessary in 
order to observe the relativistic electrons in the entire downstream region. 

In our two-dimensional purely perpendicular shock simulations with the out-of-plane 
magnetic field configuration, only drift waves along the shock surface is allowed to exist. 
Growth of the drift waves (e.g., the lower hybrid drift instability) is observed in the 
super-high Mach number shock (not shown), however, relationships between the waves and 
the electron heating or acceleration were not clearly found in the present study. A similar 
point of view must be introduced if one examines the in-plane magnetic field configuration 
in which electromagnetic waves are allowed to propagate along the shock surface. In 
supercritical shocks, two-dimensional hybrid simulations have revealed ion-scale rippled 
structures at the shock surface associated with strong electromagnetic pe rturbations, and 



the r i pple is observa ble only in the in-plane magnetic field configuration (ILowe fc Burgess 



2003 



Burgess 



20061 ) . In this case, the angle between the shock normal and the magnetic 
field direction changes locally, thereby the reflection rate of the ions is enhanced in some 
locations feeding more energy to the Buneman instability. Such a conn ection between the 



20091). 



rippled shock structure and the electron acceleration has been reported (lUmeda et al. 
The discussions are, however, limited in relatively low Alfven Mach number and small mass 
ratio cases, and understanding how the ion-scale rippled structure and the electron-scale 
Buneman instability coexist in high Mach number shocks are remained as a next step. 

Note that a two-dimensional simulation can deal with either the two-dimensional 
electron surfing acceleration in the out-of-plane magnetic field configuration or the 
rippled structure in the in-plane configuration. In latter case, the surfing acceleration 
occurs one-dimensionally and the electrons cannot escape from the trapping region. The 
acceleration efficiency will be different from the present out-of-plane case. Therefore, a 
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unified understanding of electron accelerations in liigli Macli number sliocks will be made 
possible only by a three-dimensional PIC simulation, which will be reported in near future 
by virtue of a peta-scale supercomputer system. 

This work was supported by Grant-in-Aid for Research Activity Istart-up 23840047. 
Numerical computations were partly carried out on Fujitsu FXl at JEDI of Japan Aerospace 
Exploration Agency (JAXA), Fujitsu HX600 at Nagoya University, and Cray XT4 at CfCA 
of National Astronomical Observatory of Japan (NAOJ). 
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Fig. 1. — Overview of a 2D shock structure in Run C obtained at VLgiT = 6. (a) Two- 
dimensional profile of B^. (b) Ion phase space (X-vx) plot. The color code shows number 
of particles in a logarithmic scale, (c) Electron phase space (X-energy) plot. The energy 
(7 — 1) in the y axis is in a logarithmic scale, (d) Electron number density profile and (e) 
strength of the electrostatic field (rotation-free-part of the electric field) in the zoomed-up 
region surrounded by the white-dashed lines in panel (a). 
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Fig. 2. — (a) Zoomed- up spatial profile of the electrostatic potential energy. The color code 
is normalized by the upstream kinetic energy, (b) Velocity distribution function of the ions 
sampled in the region in panel (a), (c) Fourier power spectrum of the electrostatic potential 
in two dimensions. Wave number in each axis {k^ and ky) is normalized by AVo/upe, where 
AVa = 0.4c 



-26- 




Fig. 3. — Energy spectra of the electron in the foot (gray) and downstream (black) regions 
obtained in Run C. Dashed line is a fitted Maxwell distribution to the downstream spectrum. 
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Fig. 4. — Trajectories of accelerated (red) and non-accelerated (cyan) electrons superposed 
on the two-dimensional profile of the electrostatic field strength. Open square locates the 
position of the particle at each time step following the trajectory back for 15 ^^^e^- A number 
located beside the open square denotes 7 — 1. Snapshots are taken at times of (a) fi^eT = 
513.5, (b) ^giT = 522.5, and (c) ^IgiT = 532.5, respectively. This figure is also available as 




Fig. 5. — Time histories of (a) the energy (7 — 1) and (b) the magnetic moment. Sohd 
and dashed hnes are for the accelerated and non-accelerated particles, respectively. Stars in 
panel (a) indicate the times in Figure HI 
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Fig. 6. — (a) Trajectory of the accelerated particle in the same format as Figure IH Snapshot 
is taken at QgeT = 580. This figure is also available as an mpeg animation in the electronic 
edition of the Astrophysical Journal. Time histories of (b) 7 — 1, (c) the magnetic moment, 
and (d) at the particle's position. Running average for 15 is applied to the magnetic 
field profile (d). 
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X (c/Op.) 

Fig. 7. — Comparisons of (a) B^, (b) Ey, and (c) — cV x Ej^ around the overshoot between 
Run A (dashed hne) and Run C (sohd hne) . The spatial profile across the shock is averaged 
in the y direction with respect to the overshoot position at each y position. 
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Fig. 8. — Comparison of the electrostatic field strength among simulation runs at ^giT = 6. 
Two-dimensional spatial profiles are shown for (a) Run A, (b) Run B, (c) Run C, and (d) 
Run D. 
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Fig. 9. — Energy spectra of the electron in the downstream region. The results are compared 
for (a) Run A (green), Run B (magenta), and Run C (black), and (b) Run C (black) and 
Run D (cyan). Dashed lines are fitted Maxwell distributions. 
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Fig. 10. — Upstream conditions of Runs A-D and the recent two-dimensional PIC simulations 
as a function of the mass ratio and the Alfven Mach number. Solid line indicates the marginal 
condition of eq. (|8]). Dashed lines indicate the marginal condition of the Buneman instability 
(eq. ([5])) for selected /3e conditions (/3e = 0.5, 4.5, 25). 
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Table 1: Upstream conditions of simulation runs and resultant electron accelerations 





Ma 


M/m 






eq.© 


eq.® 


ele. accel. 


Run A 


14.4 


25 


0.5 


10 


2.1 


10.3 


BI 


Run B 


16.2 


100 


0.5 


10 


4.2 


25.8 


weak 


Run C 


30.0 


100 


0.5 


10 


4.2 


25.8 


BI + adiabatic 


Run D 


30.0 


100 


4.5 


10 


12.7 


25.8 


weak 



